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1.0 


INTRODUCTION 


Recently,  Increased  Interest  has  been  directed  toward  obtaining 
better  understanding  of  the  nature  of  electron  density  irregularities  in  the 
F-reglon  of  the  high  latitude  ionosphere.  Two  areas  of  particular  Interest 
are  the  mechanisms  by  which  these  irregularities  are  created  and  the  Doppler 
spectral  characteristics  that  these  irregularities  introduce  when  scattering 
radiowaves  in  the  HF  (3-30  MHz)  frequency  band.  It  is  felt  that  these 
characteristics  may  be  useful  in  understanding  the  mechanisms  that  produce  the 
irregularities.  Also,  knowledge  of  the  Doppler  spectral  characteristics  is 
essential  for  assessing  the  clutter  effects  that  these  radars  might  have  on  HF 
radar  systems. 

During  the  past  year,  several  Doppler  spectral  studies  have  been 
conducted  using  the  Air  Force  bistatlc  research  radar  at  Ava  and  Verona,  New 
York;  a  French  scientific  radar  at  Lycksele,  Sweden;  and  the  NSF/NOAA  HF 
sounder  at  Fairbanks,  Alaska.  All  of  these  facilities  utilize  different 
Doppler  spectral  detection  schemes,  and  each  have  advantages  and  disadvan¬ 
tages.  In  this  report  we  describe  the  different  techniques  and  results  from 
each  of  these  facilities.  We  also  describe  a  new  Instrument  utilizing  a 
multipulse  technique  that  should  become  operational  within  a  few  months.  The 
multipulse  techniques  will  be  Implemented  on  a  radar  operating  in  a  selectable 
fixed  frequency  pulse  mode;  however,  it  may  equally  well  be  implemented  on  a 
swept-frequency  radar  such  as  the  Ava-Verona  facility  through  some  modifi¬ 
cation  of  that  system. 

During  the  past  year  an  appreciable  effort  was  also  devoted  toward 
evaluating  a  new  technique  for  decreasing  the  detrimental  effects  that  broad- 
banded  clutter  spectra  have  on  HF  radar  systems.  The  elements  of  this 
technique  are  described  in  the  1980  Final  Report  entitled  Doppler  Spectral 
Characteristics  of  High  Latitude  Ionospheric  Irregularities  (Greenwald, 
1980).  Results  of  this  effort  are  described  in  this  report. 
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2.0  DOPPLER  TECHNIQUES  USED  IN  STUDYING  HIGH  LATITUDE  F-REGION  IRREG- 

ULARITIES 


In  this  report  we  consider  four  different  techniques  for  studying 
the  spatial  and  temporal  dependence  of  Doppler-spectral  characteristics 
associated  with  HF  radar  backscatter  from  high  latitude  F-region 
irregularities.  The  techniques  are: 

1)  FM-CW 

2)  repetitive  pulse 

3)  variable-lag  double  pulse 

4)  multiple  pulse 

The  FM-CU  technique  is  the  method  presently  used  by  the  Ava-Verona 
blstatic  radar.  It  utilizes  repetitive  synchronized  linear  frequency  sweeps 
(quadratic  phase  variation)  in  both  the  transmitter  and  receiver.  If  the 
transmitter  and  receiver  sweep  together,  then  a  signal  transmitted  at 
frequency  F  and  scattered  from  range  Rq  would  return  to  the  receiver  when  it 
is  tuned  to  a  central  frequency  of  Fq  +  2BRq/C’.  Here  B  is  the  sweep  rate  of 
the  transmitter  and  receiver  in  units  of  megahertz  per  second  and  C’  is  the 
average  group  velocity  of  the  radio  frequency  signal  along  the  path.  The  fre¬ 
quency  offset  associated  with  the  distance  to  the  scattering  element  is  inde¬ 
pendent  of  the  instantaneous  operating  frequency;  thus  the  demodulated  signal 
at  the  output  of  the  receiver  has  a  frequency  spectrum  determined  by  the  rela¬ 
tive  strength  of  the  backscattered  signal  as  a  function  of  range. 

With  the  FM-Oi/  method  one  uses  Fourier  analysis  to  obtain  the 
scattered  signal  level  as  a  function  of  range.  A  detail  description  of  this 
process  and  the  ambiguity  function  associated  with  it  is  given  in  Rihaczek 
(1969).  For  the  purpose  of  this  report  it  is  sufficient  to  note  that  the 
complex  Fourier  coefficients  for  each  frequency  band  are  nearly  equivalent  to 
the  range  dependent  quadrature  outputs  of  a  phase  coherent  receiver  in  a 
pulsed  radar.  The  single  difference  is  that  scatterers  moving  with 

appreciable  velocity  are  Doppler  shifted  to  appear  as  if  they  were  located  at 
a  different  range. 
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The  analogy  between  repetitively  swept  FM-CW  radars  and  repetitively 
pulsed  fixed  frequency  radars  may  be  carried  one  step  further.  Let  us  assume 
that  there  is  a  stationary  scatterer  such  as  a  mountain  at  range  Rq.  Then  at 
a  group  delay  x  •  2Rq/C',  the  quadrature  outputs  of  a  phase  coherent  receiver 
would  maintain  constant  values.  In  a  like  manner,  If  the  mountain  were  to 
scatter  a  signal  from  an  FM-CW  radar,  then  the  complex  Fourier  coefficient 
corresponding  to  the  frequency  2BRg/C'  and  range  Rq  would  also  remain  constant 
with  time.  Alternatively,  If  the  scatterer  at  range  Rq  were  moving,  then  both 
the  quadrature  receiver  outputs  at  group  delay  x  and  the  complex  Fourier 
coefficient  at  frequency  2BRq/C'  would  be  changing  with  time  and  they  would 
change  In  an  identical  manner. 

We  can  now  consider  a  series  of  consecutive  pulses  from  a  fixed  fre¬ 
quency  radar  or  consecutive  sweeps  of  an  FM-CW  radar.  From  these  one  can 
obtain  a  series  of  consecutive  samples  of  the  output  of  the  phase  coherent 
receiver  at  group  delay  x  or  a  series  of  consecutive  values  for  the  complex 
Fourier  coefficient  corresponding  to  frequency  2BRq/C'.  Either  of  these 
series  may  be  taken  as  a  time  series  which  may  be  Fourier  transformed  to  yield 
the  Doppler  spectrum  of  the  scatterers  located  at  range  Rq. 

When  performing  Doppler  spectral  analysis,  certain  restrictions  must 
be  taken  into  account.  Firstly,  the  Doppler  width  and  mean  Doppler  velocity 
of  the  scatterers  at  range  Rq  determine  the  repetition  frequency  of  the  FM-CW 
sweeps  and  the  fixed  frequency  radar  pulses.  Ideally  these  should  be  In  ex¬ 
cess  of  twice  the  maximum  frequency  in  the  scatterer' s  Doppler  spectrum  for  a 
given  radar  wavelength.  For  example,  if  the  scatterers  possess  velocities  of 
1000  m/s  relative  to  a  radar  operating  at  15  MHz,  then  the  repetition  frequen¬ 
cy  should  be  in  excess  of  200  Hz. 

Having  selected  a  repetition  frequency  for  the  radar  pulse  or  sweep 
rate,  one  is  immediately  confronted  with  the  second  restriction.  With  high 
repetition  frequencies,  it  Is  entirely  possible  that  the  receiver  will  receive 
backscatter  from  multiple  ranges  simultaneously.  For  a  fixed  frequency  radar 
with  a  repetitive  frequency  of  200  Hz,  scatterers  at  a  group  delay  x  would  be 
superimposed  upon  scatterers  at  group  delays  of  x+  5  ms,  x  +  10  ms,  etc.  In  a 
like  manner  for  a  FM-CW  system,  scatterers  located  at  Rq/ C  would  be 
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superimposed  upon  scacterers  at  Rq/ C*  +  5  ms,  etc.  The  severity  of  this 
problem  Is  ultimately  dependent  on  three  parameters:  1)  the  radar  operating 
frequency,  2)  the  radial  velocity  distribution  of  the  scatterers,  and  3)  the 
range  Interval  over  which  backscatter  Is  observed. 

As  we  have  described  Doppler  spectral  studies  using  repetitively 
pulsed  radars  along  with  our  description  of  FM-CM  radars,  there,  is  no  need  to 
describe  this  approach  further.  It  should,  however,  be  noted  that  the  French 
SAFARI  radar  located  In  Lycksele,  Sweden  operates  in  a  repetitively  pulsed 
selectable  fixed-frequency  mode. 

The  variable  lag  double-pulse  technique  is  an  approach  in  which  one 
may  obtain  unallased  measurements  of  large  Doppler  velocities  while  simultan¬ 
eously  avoiding  the  problem  of  superimposing  backscatter  from  a  number  of  dif¬ 
ferent  ranges.  This  approach  Involves  the  determination  of  the  complex  auto¬ 
correlation  function  rather  than  the  Doppler  spectrum.  However,  once  the 
autocorrelation  function  Is  obtained  the  Doppler  spectrum  may  be  determined  by 
Fourier  transformation. 

We  can  understand  the  double  pulse  approach  as  follows :  consider 
two  pulses  emitted  by  a  transmitter  and  separated  by  time  t.  They  are  scat¬ 
tered  from  a  target  at  range  Rq.  At  time  delay  t  after  each  pulse  a  back- 
scattered  response  will  appear  at  the  receiver.  Initially,  let  us  assume  the 
target  to  be  discrete  and  moving  toward  or  away  from  the  radar  with  constant 
velocity.  The  sampled  outputs  of  the  quadrature  receiver  may  be  denoted  as: 


C^(t)  -  A^(t)  +  IBj^(t) 

2.1 

C2(T,t)  -  AjCr.t)  +  IB^Cx.t) 

2.2 

where  the  subscripts  1  and  2  denote  the  sample  associated  with  the  first  and 
second  transmitter  pulse  at  group  delay  t.  The  autocorrelation  func¬ 
tion  R(T,t)  Is  given  by 
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It  possesses  real  and  Imaginary  parts  given  by 


X(T,t)  =  Aj^(T)A2(T,t)  +  Bj(T)B2(T,t)  “  a  co6(2vt/A)  2.4 

and 

Y(T,t)  -  A2(T)B^(T,t)  -  Aj(T)B2(T,t)  -  a  sln(2vt/X)  2.5 


where  a  Is  a  constant,  v  Is  the  radial  velocity  of  the  target  and  X  Is  the 
wavelength  corresponding  to  the  radar  operating  frequency.  One  can  see  that 
the  real  and  Imaginary  parts  of  the  autocorrelation  function  trace  out  sinu¬ 
soids  that  are  In  quadrature.  If  these  quantities  are  determined  at  discrete 
lag  Increments,  t=ntQ,  between  the  two  transmitted  pulses,  then  the  resulting 
values  could  be  Fourier  transformed  to  yield  a  line  spectrum  at  f  =  ±2v/X  . 
The  sign  of  the  Doppler  shift  Is  determined  by  whether  the  target  is  moving 
toward  or  away  from  the  radar. 

In  the  case  of  a  distributed  target  or  one  moving  with  variable 
velocity  the  analysis  is  slightly  more  complex.  The  problem  is  described  In 
greater  detail  by  Rummler  (1968).  For  the  present  report  it  Is  sufficient  to 
note  that  the  autocorrelation  functions  will  generally  appear  as  damped 
sinusoidal  functions. 

We  can  now  consider  some  of  the  advantages  of  the  variable-lag 
double  pulse.  Since  we  have  only  transmitted  a  double  pulse  rather  than  a 
repetitive  string  of  pulses,  at  any  time  the  receiver  will  receive  the  back- 
scattered  return  from  no  more  than  two  ranges.  Moreover,  In  the  cross  product 
calculation  only  one  range  will  contribute  to  an  average  non-zero  result. 
This  can  be  seen  by  allowing  2. 1-2. 2  to  be  generalized  so  that  there  are  tar¬ 
gets  at  different  ranges.  Then 


S  "  Ai(T2.t)  +  IBiCr^.t)  +  A2(Ti,t)  +  iB^Crpt) 


where  the  quadrature  receiver  output  Ci  contains  backscattered  contributions 
from  the  first  transmitted  pulse  at  group  delay  and  from  the  second  trans¬ 
mitted  pulse  at  group  delay  Tq<  Ti  ,  and  C2  contains  contributions  from  the 
first  transmitted  pulse  at  group  delay  12  >  and  from  the  second  transmitted 
pulse  at  group  delay  ti.  Under  the  generally  valid  assumption  that  back- 
scatter  from  distinct  ranges  is  statistically  Independent,  one  can  see  that 


<R(Ti,t)>  -  <CiC2> 


”  <[Ai(ti)  +  1Bi(ti)]  [A2(Ti,t)  +  lB2(Ti,t)]> 


where  <  >  denotes  the  expectation  value.  Note  that  this  result  is  essentially 
independent  of  the  spacing  between  the  two  transmitted  pulses.  Hence,  very 
large  Doppler  velocities  may  be  correctly  analyzed  without  superimposing 
contributions  from  distinct  ranges.  In  this  operational  mode  one  must  only 
take  care  that  the  spacing  between  double  pulse  pairs  is  sufficiently  large  so 
backscatter  is  received  from  only  one  pulse  pair  at  a  time. 


In  terms  of  comparison  it  is  possible  to  consider  the  equivalent 
autocorrelation  function  that  results  from  a  repetitive  transmitter  pulse 
sequence.  One  can,  for  example,  rewrite  (2. 6-2. 7)  as 


C,  =  A„(tJ  +  iB„  (tJ  +  A,(t,)  +  iB.  (t.)  +  A„(t„)  +  1B„(t.)  2.9 


Aj(T2,t)  +  lBj(T2ft)  +  A^CTj.t)  +  iB2(Tj,t)  +  A^Cr^.t)  +  IB^Ct^.t) 


where  subscripts  0-3  on  the  A's  and  B's  denote  four  sequential  pulses  in  a 
repetitive  pulse  series.  The  associated  autocorrelation  function  is  given  by 


<R(T,t)>  -  <C^C*> 

=  <(Aq(t2)  +  1Bq(t2)I 

+  <[Aj(t^)  +  IBj(Tj^)] 
+  <IA2(Tq)  +  1B2(Tq)] 


lA^(T2,t)  +  lBj(T2,t))> 

[A2(Ti,t)  +  iB2(Tj,t)]> 

lAgCr^.t)  +  lB2(TQ,t)l>. 


2.10 


Thus,  by  using  repetitive  pulse  we  have  superimposed  the  results  from  three 
distinct  group  delays,  T2.  and  Tq  . 

In  fairness  it  should  be  noted  that  the  variable  lag  double  pulse  is 
rather  inefficient  in  terms  of  maintaining  a  high  average  power  output.  This 
is  especially  true  for  the  longer  lags  of  the  autocorrelation  function  for 
which  the  transmitter  is  generally  inactive.  The  approach  also  assumes  that 
the  scattering  process  is  time  stationary  during  the  Interval  that  the  various 
lags  of  the  autocorrelation  function  are  being  calculated.  Often  and  espe¬ 
cially  in  the  high  latitude  ionosphere  this  is  not  the  case. 

The  variable  lag  double  pulse  approach  has  been  used  for  spectral 
measurements  with  the  STARE  radar  system  in  Scandanavla  and  with  the  NSF/NOAA 
HF  sounder  at  Fairbanks,  Alaska.  With  the  latter  instrument  difficulties  due 
to  the  non-stationarity  of  the  scattering  process  were  readily  apparent. 

A  vast  Improvement  over  the  simple  variable  lag  double  pulse  scheme 
is  the  multipulse  approach.  In  this  technique  one  transmits  a  number  of 
pulses  with  carefully  selected  spacings.  The  criteria  for  choosing  the  pulse 
pattern  are:  1)  that  one  should  obtain  with  one  pulse  burst  as  many  lags  of 
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an  autocorrelation  function  as  possible  and  2)  that  one  should  avoid  pulse 

schemes  that  allow  superposition  of  different  group  delays.  Farley  (1972)  has 
presented  a  number  of  possible  schemes  for  Incoherent  scatter  measurements. 
For  HF  measurements  In  the  high  latitude  ionosphere  the  problem  Is  somewhat 
less  restrictive,  and  It  Is  possible  to  choose  pulse  schemes  that  Farley  may 

not  have  considered.  In  particular,  the  pulse  scheme  shown  in  Figure  2.1  is 

one  that  is  presently  being  implemented  on  a  HF  radar  for  high  latitude 
irregularity  studies. 

The  pulse  pattern  shown  in  Figure  2.1  utilizes  seven  transmitted 
pulse  with  spacings  of  2tQ,  tQ,  ^*’0  *"0*  This  combination 

contains  all  possible  lags  for  zero  to  Ibtg.  If  one  examines  the  possible 
lags  from  this  sequence  one  sees  that  there  are  two  possible  group  delays  that 
will  contribute  to  lags  of  tQ,  2tQ  and  3tQ.  However,  these  group  delays  are 

separated  by  lOtQ.  Assuming  a  minimum  value  of  Iq  =  2ms,  one  sees  that  the 

range  ambiguity  occurs  only  if  there  is  concurrent  scatter  from  regions  with 

group  delays  differing  by  at  least  20  ms.  This  is  an  extremely  unlikely  event 

even  for  HF  backscatter  and  thus  it  may  be  safely  ignored.  There  also  is  a 
range  ambiguity  for  a  lag  of  13tQ  with  three  distinct  group  delays  contri¬ 
buting  to  the  cross  product.  However,  as  there  should  be  relatively  little 
residual  correlation  remaining  at  this  long  lag,  the  potential  confusion  may, 
with  caution,  be  ignored. 

With  the  multipulse  technique  there  is  a  considerable  increase  in 
average  transmitted  power.  If  we  consider  a  sixteen  lag  autocorrelation 
function  then  the  time  required  to  transmit  a  complete  16  lag  multipulse 
pattern  is  16tg.  In  contrast,  if  one  uses  variable  lag  double  pulse  methods, 
the  complete  transmission  would  require  128tQ.  In  addition,  for  a  multipulse 
pattern  one  needs  to  allow  for  the  group  delay  to  the  scattering  region  only 
once,  whereas  for  the  double  pulse  approach  one  must  include  this  group  delay 
for  each  lag  of  the  correlation  function.  Thus  the  multipulse  approach  is  at 
least  eight  times  as  efficient  as  the  variable  lag  double  pulse  approach  and 
its  efficiency  increases  approximately  linearly  with  the  number  of  lags  in  the 
correlation  function. 
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A  shortcoming  of  the  multipulse  approach  Is  that  there  Is  Increased 
noise  due  to  uncorrelated  scatter  from  unwanted  ranges  (e.g.  contribution  such 
as  those  In  Equation  2. 6-2. 7).  Fortunately  for  HF  studies,  the  entire  trans¬ 
mitted  pulse  generally  does  not  cause  concurrent  backscatter  at  the  receiver 
and,  therefore,  the  problem  Is  considerably  less  severe  that  It  might  appear. 

There  Is  one  final  and  Important  point  with  respect  to  the  spectral 
and  correlation  techniques  described  In  this  section.  Specifically,  It  Is 
entirely  possible  to  combine  the  FM  or  quadratic  phase  techniques  with  the 
double  pulse  or  multiple  pulse  techniques.  Under  these  circumstances  a  pulse 
would  be  replaced  with  a  linear  FM  sweep  with  duration  less  than  tQ.  Also, 
control  must  be  Included  to  determine  the  occurrence  and  spacing  of  the  FM 
sweeps.  One  would  then  calculate  a  Fourier  transform  for  each  sweep  to  obtain 
the  complex  Fourier  coefficients  corresponding  to  each  group  delay  and  use 
these  coefficients  In  (2.8)  In  the  same  manner  as  one  would  use  the  quadrature 
outputs  of  a  phase  coherent  receiver.  This  combination  would  produce  a  rather 
efficient  transmission  made  with  relatively  high  averge  power  and  no  range 
ambiguities  or  Doppler  spectral  aliasing. 


3.0 


DOPPLER  SPECTRA  FROM  HIGH  LATITUDE  P-REGION  IRREGULARITIES 


Having  considered  a  variety  of  methods  for  the  determination  of 
Doppler  spectra  and  autocorrelation  functions  from  high  latitude  F-region 
Irregularities,  we  now  examine  some  typical  results.  Figures  3. 1-3. 6  repre¬ 
sent  examples  of  Doppler  spectra  observed  with  the  Ava-Verona  FM-CW  system 
during  the  evening  hours.  The  measurements  were  made  on  27  January  1981 
between  0100-0300  UT  (2000-2200  local  time).  For  these  measurements  the  radar 
was  swept  over  the  frequency  range  16.387-16.389  Mlz  with  a  sweep  rate  of  100 
KHz/s  and  the  receiving  antenna  was  directed  toward  geographic  north  (approx¬ 
imately  10“  to  the  west  of  geomagnetic  north) .  The  times  given  represent 
universal  time.  The  range  Index  represents  the  range  interval  for  which  the 
Doppler  spectrum  is  analyzed.  For  these  examples  the  range  is  determined  by 
multiplying  the  index  value  by  75  km.  In  the  case  of  Figure  3  the  range  cor¬ 
responds  to  1350  km. 

Given  the  operational  parameters,  the  Doppler  spectra  in  Figures 
3. 1-3. 6  have  an  aliasing  frequency  of  25  Hz  which  corresponds  to  a  radial 
velocity  of  approximately  230  m/s.  Ignoring  for  the  moment  the  narrow  spec¬ 
tral  peak  at  +25  Hz  (this  will  be  discussed  in  the  next  section),  one  sees 

that  the  Doppler  spectra  are  generally  quite  broad.  The  narrowest  Doppler 
width  is  observed  in  Figure  3.2  where  it  is  aproxlmately  160  m/s  and  the 
widest  is  observed  in  Figure  3.6  where  the  spectrum  is  virtually  flat. 

Clearly,  most  of  the  spectra  are  rather  strongly  frequency  aliased  and  have 
spectra  widths  in  excess  of  250  m/s. 

One  can  attempt  to  locate  the  peak  of  each  of  the  spectra;  however, 
again  due  to  the  frequency  aliasing,  the  values  may  not  be  correct  even  with 
regard  to  the  sign  of  the  shift.  This  is  due  to  the  fact  that  spectral 

aliasing  beyond  the  Nyqulst  frequency  comes  in  at  the  negative  Nyquist  fre¬ 
quency  as  can  be  seen  in  Figure  3.5.  Thus  the  evidence  in  the  data  of  sign 
changes  in  the  mean  Doppler  shift  is  quite  possibly  due  to  changes  in  the 
magnitude  of  a  unidirectional  drift  velocity. 

In  an  attempt  to  obtain  unallased  spectra  an  experiment  was  con¬ 

ducted  on  10  April  1981  with  a  sweep  repetition  frequency  of  200  Hz  and  an 
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average  operating  frequency  of  18.5  Hlz.  For  this  study  a  Doppler  velocity  of 
810  m/s  would  have  produced  frequency  aliasing.  Unfortunately,  the  resulting 
Doppler  spectra  were  quite  similar  to  Figure  3.1,  ^rf^lch  again  suggested  fre¬ 
quency  aliasing.  No  further  tests  were  conducted  at  these  high  sweep  repeti¬ 
tion  frequencies. 

It  Is  worthwhile  comparing  these  Doppler  observations  with  some  the¬ 
oretical  expectations  as  to  the  motions  of  F- region  Irregularities.  Greenwald 
(1980)  has  postulated  that  F-reglon  Irregularities  might  drift  under  the  in¬ 
fluence  of  the  Ionospheric  electric  field.  In  this  case,  the  drift  at  these 
local  times  would  be  westward  with  a  magnitude  of  a  few  hundred  to  more  than 
one  thousand  meters  per  second.  As  the  radar  Is  directed  to  the  west  of  geo¬ 
magnetic  north,  a  negative  Doppler  shift  Is  expected.  In  the  figures  positive 
Doppler  shifts  are  at  times  observed.  While  this  may  indicate  an  inconsistency 
with  the  electric  field  drift  hypothesis,  it  is  also  possible  that  it  is  in¬ 
dicative  of  spectral  aliasing.  Clearly,  better  spectral  results  are  required. 

The  French-Swedish  SAFARI  radar  (Hanulse  et  al.,  1981)  has  now  been 
operating  for  one  year.  In  this  period,  a  number  of  campaigns  have  been  con¬ 
ducted  and  results  similar  to  those  shown  In  Figure  3.7  have  been  obtained.  In 
this  example,  the  radar  was  operated  In  the  14  MHz  amateur  band  with  a  5”  hor¬ 
izontal  beamwldth  directed  toward  geographic  north  (~  12°  to  the  east  of 
corrected  geomagnetic  north).  The  radar  was  operated  In  a  pulse  mode  with  a 
pulse  repetition  frequency  of  LOO  Hz. 

One  can  see  In  Figure  3.7  that  with  a  pulse  repetition  frequency  of 
100  Hz  there  was  relatively  little  frequency  aliasing  of  the  spectra.  The 
spectral  widths  varied  from  8-24  Hz  (80-250  m/s)  and  the  mean  Doppler  shifts 
ranged  from  200-380  m/s.  At  the  local  time  of  this  measurement  (~  2315)  the 
electric  field  Induced  drift  would  be  In  a  southwestward  direction.  The  sign 
of  the  Doppler  shift  Is  consistent  with  this  expected  drift;  however,  with 
only  a  single  Doppler  component  It  Is  difficult  to  estimate  what  the  total  Ir¬ 
regularity  drift  speed  is.  It  should  be  noted  that  there  are  range  dependent 
variations  In  the  observed  Doppler  velocity  which  may  be  associated  with  spa¬ 
tial  variation  of  magnitude  and/or  direction  of  an  electric  field  Induced 
drift. 
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ELIMINATION  OF  SPREAD  SPECTRA 


In  a  recent  report  Greenwald  (1980)  suggested  that  cross  spectral 
analysis  might  be  used  to  eliminate  spread  Doppler  spectra  associated  with 
Ionospheric  Irregularities  while  retaining  discrete  spectra  that  might  be 
associated  with  aircraft.  The  essence  of  this  approach  was  to  calculate  the 
time-averaged  cross  spectrum  of  two  lagged  data  sets  for  which  the  lag  time 
was  greater  than  the  decorrelation  time  of  the  ionospheric  process  contri¬ 
buting  to  the  broad  spectral  component. 

Using  Equations  5.9  -  5.16  from  Greenwald  (1980)  we  can  write 

S'(R,aj)  -  A(o))  B*(u))  4.1 

as  the  cross  spectrum  of  the  signal  that  Is  backscattered  from  range  R.  One 
may  express  the  complex  Fourier  coefficients  A  (u)  and  B*(m)  in  terms  of  dis¬ 
cretely  sampled  time  series.  These  are  given  by 


A(R,n,p)  X(R,j,ji)W(j)exp(-2nljn/N)  4.2 

B(R,n,p)  jiJ  X(j+C,p)W(j)exp(-27rljn/N)  4.3 

where  X(R,j,vi)  is  the  complex  demodulated  output  of  the  receiver  at  the  group 

delay  corresponding  to  range  R,  W(j)  Is  the  data  window  (Blackman  and  Tukey, 

1939),  and  y  Is  the  number  of  the  data  set  when  multiple  data  sets  are  being 

transformed.  Assuming  N  samples  per  data  set  and  a  sampling  frequency  f„, 

0 

each  data  set  Is  obtained  In  the  time 

AT  -  N/f  4.4 

s 

while  the  lag  time  between  the  data  sets  used  in  the  cross  spectral  determi¬ 
nation  Is 


T  -  C/fg 


4.5 
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The  cross  spectral  estimator  Is  then  given  by 


Pjjjj  (n)  -  A(n,p)B*(n,u)  4.6 

and  the  variance  of  this  estimator  is 

o(n,M)  -  S(R,n)//M  4.7 

where  M  is  the  number  of  spectra  that  are  averaged  and  S(R,n)  Is  the  spectral 
power  in  frequency  Interval  n  of  the  signal  backscattered  from  range  R.  It 
has  a  spectral  estimator  given  by 


A(n,p)A*(n,g) 


4.8 


If  the  lag  time  T  is  sufficiently  large  that  A(n,  y)  and  B(n,u)  are 
uncorrelated,  then  for  M  +  *,  the  cross  spectrum  will  average  to  zero  incoher¬ 
ently.  That  is 


Pxx  (n)  =  ~  a(n,M)  4.9 

c 

In  contrast  contributions  to  A(n,u)  and  B(n,y)  that  remain  coherent  after  lag 
time  T  will  maintain  the  same  power  level  In  the  cross  spectrum  as  in  the 
power  spectrum. 

While  the  results  to  be  presented  in  this  section  support  these  con¬ 
clusions,  they  also  bring  out  a  point  that  had  been  overlooked  in  the  previous 
report  (Greenwald,  1980).  Specifically,  these  results  show  that  ones  ability 
to  discern  a  discrete  target  within  a  spread  spectrum  is  controlled  by  the 
variance  of  the  spectral  estimators.  Since  the  variance  of  both  the  power 
spectrum  and  cross  spectrum  decrease  at  a  rate  that  is  approximated  by  (4.7), 
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one  method  does  not  have  any  significant  advantage  over  the  other.  In  prac¬ 
tice,  the  cross  spectral  approach  does  appear  to  have  a  variance  that  is  3  dB 
less  than  the  power  spectral  approach.  This  difference  may  be  due  to  the  fact 
that  both  the  in-phase  and  quadrature  portions  of  the  cross-spectrum  are  being 
averaged  in  cross-spectral  calculation. 

During  the  past  year,  several  studies  have  been  conducted  to  eval¬ 
uate  the  usefulness  of  the  cross  spectral  approach.  The  narrow  spectral  peak 
located  at  the  right  hand  edge  of  the  Doppler  spectra  in  Figures  3.1  -  3.6 
represent  an  additional  discrete  frequency  Inserted  into  the  receiver  to 
simulate  the  presence  of  a  discrete  target.  These  figures  represent  the 
average  of  100,  200  and  400  Doppler  power  spectra  and  one  sees  that,  aside 
from  the  added  signal,  the  variance  in  the  spectral  estimates  is  typically  of 
the  order  of  that  which  one  would  expect. 

The  same  data  have  been  used  to  calculate  the  cross  spectral  power 
within  the  Doppler  spectra  for  a  lag  time  of  1.28  s.  That  is,  instead  of  cal¬ 
culating  the  Doppler  power  spectrum,  two  sets  of  Fourier  coefficients  obtained 
1.28  s  apart  were  cross  multiplied  to  yield  the  Doppler  cross  power  spec¬ 
trum.  The  results  are  shown  in  Figures  4. 1-4. 6.  In  these  examples,  one  can 
see  the  reduction  of  the  overall  level  of  the  Doppler  spectrum.  In  Figures 
4.1  and  4.6  it  is  roughly  10%  of  the  levels  in  Figure  3.1  and  3.6.  In  Figures 
4.2-4 .4  for  which  200  spectra  were  averaged,  the  level  is  approximately  7%  of 
the  level  in  Figures  3.2-3 .4.  Finally  in  Figure  4.5  for  which  400  spectra 
were  averaged,  the  level  of  the  background  spectrum  is  only  5%  of  that  in 
Figure  3.5  (note  the  change  in  the  vertical  axis  scaling  for  Figure  4.5). 

If  one  compares  the  discrete,  inserted  frequency  in  Figures  3. 1-3. 6 
and  4. 1-4. 6,  one  sees  that  there  is  no  measureable  reduction  in  this  component 
as  a  result  of  cross  spectral  analysis.  Any  reduction  in  the  frequency  band 
containing  the  Inserted  frequency  is  a  reduction  in  the  ionospheric  signal  on 
which  the  inserted  frequency  is  sitting. 

Again,  comparing  Figure  3. 1-3. 6  and  4. 1-4. 6,  one  sees  that  in  all 
cases  the  discrete  spectrum  is  easier  to  Identify  in  the  latter  spectra  than 
in  the  former.  Also,  the  variance  of  the  spectral  estimators  of  the  cross 


spectra  appears  Co  be  less  chan  chose  of  Che  power  spectra  as  has  been  noted 
previously.  However,  In  all  cases  for  which  the  discrete  frequency  Is  readily 
Identifiable  In  the  cross  spectrum,  I  believe  It  could  also  be  Identified  In 
the  power  spectrum  albeit  with  greater  difficulty. 

In  passing  it  should  be  noted  that  cross  spectral  analysis  may  also 
be  used  for  the  first  Fourier  transform  of  the  FM-CW  technique.  When  used  for 
this  purpose  It  Is  also  capable  of  discerning  discrete  targets.  Figure  4.7 
Illustrates  the  range  dependence  of  the  backscaCtered  signal  on  Day  346  at 
0344  UT.  In  this  example,  each  100  Hz  corresponds  to  a  range  Interval  of  150 
km  with  zero  Hertz  having  been  offset  to  a  range  of  1200  km  and  the  spectra 
have  been  plotted  on  a  log  scale.  The  spectral  peaks  observed  at  Integral 
multiples  of  50  Hz  represents  a  modulation  of  the  backscaCtered  signal  by  Che 
sweep  repetition  frequency  of  the  transmitter. 

The  cross  power  spectrum  shown  In  Figure  4.8  utilizes  the  same  data 
as  has  been  used  to  obtain  Figure  4.7;  however,  the  spectra  have  been  calcu¬ 
lated  with  a  1  second  lag.  Here  one  sees  that  the  signal  level  Is  approxi¬ 
mately  15  dB  less  than  that  In  Figure  4.7  commensurate  with  the  nearly  1000 
spectra  Chat  have  been  averaged  Co  obtain  these  results. 

One  final  point  should  be  made  with  regard  to  both  power  or  cross 
spectral  analyses.  If  a  transmitter  has  a  pulse  or  sweep  repetition  frequency 
of  64  Hz  and  one  Is  calculating  64  sample  Doppler  spectra,  then  one  second  Is 
required  to  obtain  the  data  for  one  spectrum.  In  order  to  average  100 
spectra,  100  seconds  of  data  Is  required.  Due  to  a  number  of  reasons  such  as 
changing  Ionospheric  conditions,  100  seconds  Is  probably  the  practical  upper 
limit  for  averaging  spectra.  The  variance  of  the  spectral  estimators  when  100 
spectra  are  averaged  is  10%  Implying  that  it  Is  unlikely  that  one  could 
Identify  a  discrete  signal  that  Is  more  than  10  dB  below  the  total  Doppler 
power  In  any  given  frequency  band. 


5.0 


SUMMARY 


In  this  report  ue  have  considered  several  techniques  used  In  FM~CW 
and  pulsed  radars  for  obtaining  the  Doppler  sectrum  associated  with  high 
latitude  F-reglon  electron  density  Irregularities  as  observed  at  HF  frequen¬ 
cies.  It  Is  found  that  one  Is  typically  confronted  with  the  conflicting 
problems  of  pulsing  or  sweeping  the  radar  at  a  sufficiently  high  repetition 
frequency  to  avoid  frequency  aliasing  of  the  Doppler  spectrum  while  at  the 
same  time  avoiding  range  aliasing  of  the  backscattered  signals.  Results 
obtained  with  the  FM-W  radar  at  Ava-Verona,  New  York  and  a  French-Swedlsh 
radar  at  Lycksele,  Sweden  Indicate  that  mean  Doppler  velocities  and  Doppler 
widths  In  excess  of  400  ms  may  not  be  uncommon.  For  radars  operating  In  the 
frequency  range  of  15-20  Mlz,  this  would  require  repetition  frequencies  of  at 
least  100  Hz.  With  these  repetition  frequencies  the  unambiguous  range 
resolution  Is  10  ms  or  1500  km.  While  there  are  periods  for  which  such 
conditions  may  be  satisfied,  at  other  times  higher  repetition  frequencies  may 
be  required  or  the  backscatter  may  be  extended  over  greater  ranges.  Under 
these  circumstances  more  sophisticated  pulse  modulation  techniques  should  be 
sought . 


One  approach  which  might  be  used  to  a  greater  degree  is  the 
multipulse  technique.  This  technique  utilizes  special  pulse  or  FM  sweep 
patterns  that  enable  one  to  determine  most  or  all  of  the  lags  of  an  auto¬ 
correlation  function.  Having  obtained  the  autocorrelation  function  one  can 
easily  determine  the  Doppler  spectrum  via  Fourier  transformation.  The 
multipulse  approach  Is  particularly  powerful  since  It  allows  one  to  obtain 
very  high  equivalent  pulse  or  sweep  repetition  frequencies  while  avoiding 
ambiguity  as  to  the  range  of  the  backscattered  signal. 

In  this  report  we  also  present  the  results  of  several  tests  on 
determining  the  usefulness  of  cross  spectral  analysis  as  a  means  of  separating 
discrete  and  spread  spectra.  While  the  results  show  that  cross  spectral 
analysis  enables  easier  identification  of  discrete  spectral  contributions 
within  a  spread  background  spectrum,  considerable  Improvements  are  not  be  be 
expected.  This  conclusion  Is  reached  by  noting  that  the  standard  deviation  of 
the  spectral  estimators  in  the  cross  spectrum  Is  only  slightly  less  than  the 
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standard  deviation  of  power  spectral  estimators.  Furthermore,  the  ratio  of 
these  standard  deviations  Is  independent  of  the  number  of  power  or  cross 
spectra  that  have  been  averaged.  Since  ones  ability  to  discern  a  discrete 
spectral  contribution  Is  ultimately  limited  by  the  Intensity  of  this  component 
with  respect  to  the  standard  deviation  of  the  spectral  estimators,  only 
nominal  Improvements  can  be  expected. 
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FIGURE  (OPTIONS 


FIGURE  2.1  Seven  transmitter  pulse  pattern  that  provides  autocorrelation 
function  lags  from  0-16  tQ.  This  pattern  has  range  ambiguity 
for  lag  13  tQ.  At  HF  frequencies  there  will,  in  general,  not  be 
any  range  ambiguity  for  lags  tQ,  2  tQ,  and  3  tQ. 

FIGURE  3.1  Doppler  spectrum  of  F-reglon  irregularities  at  a  range  of  1350 
km  at  0105  UT  on  27  January  1981  (Ava-Verona  radar) . 

FIGURE  3.2  Doppler  spectrum  of  F-reglon  irregularities  at  a  range  of  1350 
km  at  0111  UT  on  27  January  1981  (Ava-Verona  radar). 

FIGURE  3.3  Doppler  spectrum  of  F-region  irregularities  at  a  range  of  1350 
km  at  0118  UT  on  27  January  1981  (Ava-Verona  radar). 

FIGURE  3.4  Doppler  spectrum  of  F-reglon  irregularities  at  a  range  of  1350 
km  at  0129  UT  on  27  January  1981  (Ava-Verona  radar) . 

FIGURE  3.5  Doppler  spectrum  of  F-reglon  irregularities  at  a  range  of  1350 
km  at  0132  UT  on  27  January  1981  (Ava-Verona  radar). 

FIGURE  3.6  Doppler  spectrum  of  F-region  irregularities  at  a  range  of  1125 
km  at  0230  UT  on  27  January  1981  (Ava-Verona  radar) . 

FIGURE  3.7  Range  dependent  Doppler  spectra  (group  delay  5.6-8 .6  ms)  of  F- 
reglon  ionospheric  irregularities  as  observed  with  the  French- 
Swedlsh  radar  at  Lycksele,  Sweden. 

FIGURE  4.1  Doppler  cross  spectrum  obtained  from  the  data  used  in  Figure 

3.1.  Discrete  frequency  Inserted  at  25  Hz. 

FIGURE  4.2  Doppler  cross  pectrum  obtained  from  the  data  used  in  Figure 

3.2.  Discrete  frequency  Inserted  at  25  Hz. 
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FIGURE  4.7  Range  Fourier  transform  of  FM-CW  receiver  output  from  the  Ava- 
Verona  radar  on  11  December  1980.  Spectral  lines  are  due  to 
modulation  of  the  received  signal  by  the  repetition  frequency  of 
the  FM-CM  sweep. 

FIGURE  4.8  Same  range  Fourier  transform  as  In  Figure  4.7  except  the  power 
profile  Is  calculated  using  cross  spectral  analysis  with  a  1 
second  lag.  Note  the  overall  15  dB  reduction  In  the  signal 

level  due  to  lack  of  coherence  over  this  lag  time. 
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